Ideal gas in nonextensive optimal Lagrange multipliers formalism 
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Based on the prescription termed the optimal Lagrange 
multipliers formalism for extremizing the Tsallis entropy in- 
dexed by q, it is shown that key aspects of the treatment of 
the ideal gas problem are identical in both the nonextensive 
g / 1 and extensive q = 1 cases. 
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I. INTRODUCTION 

The purpose of this work is to discuss the classical 
ideal gas problem in Tsallis thermostatistics within the 
framework of the method of optimal Lagrange multipliers 
(OLM) recently proposed in Ref. [0. 

The ideal gas problem in the normalized Tsallis ther- 
mostatistics has recently been exhaustively discussed in 
Ref. Q. Correlations induced by Tsallis' nonextensivity 
were analyzed in Ref. ||. An interesting relation, analo- 
gous to that yielding the mean energy of the (ordinary) 
ideal gas was thereby found. However, the ensuing re- 
sults still depend upon Tsallis' nonextensivity index q, 
even if one deals with a purely classical system. Here, we 
show that more insights can be gained into the problem 
if i) one takes advantage of a degree of freedom implicit 
in Tsallis' formalism, i.e., choice of the entropy constant 
kx, playing the role of Boltzmann's constant fee, and ii) 
uses the OLM technique in Ref. Q . 

Tsallis' thermostatistics [4-18] is by now known to offer 
a nonextensive generalization of traditional Boltzmann- 
Gibbs statistical mechanics. A key ingredient in this for- 
malism is the introduction of a particular definition of 
expectation value termed the normalized g-expectation 
value Q. Actually, during the last ten years before the 
work in Ref. , several proposals have been made regard- 
ing definition of expectation value j^]. No matter what 
definition one chooses, ordinary Boltzmann-Gibbs results 
are always reproduced in the extensive limit q — * 1. 

Tsallis thermostatistics involves extremization of Tsal- 
lis' entropy 



by recourse to Lagrange's constrained variational tech- 
nique. Here, x is a phase space element (N particles in a 
D-dimensional space), q € 5i is Tsallis' nonextensivity in- 
dex, / stands for any normalized probability density, and 
kx is the entropy constant that is akin to the celebrated 
Boltzmann constant appearing in traditional statistical 
mechanics. 

The formalism of Ref. gives rise to a non-diagonal 
form of the Hessian associated with examining the ex- 
tremum structure of the entropy and free energy through 
the Legendre transform procedure. In Ref. a new ap- 
proach has been advanced, in which the Hessian is diag- 
onal. This approach was shown to enormously facilitate 
ascertaining what kind of extrema the Lagrange tech- 
nique yields. 

A quite interesting fact is to be emphasized concern- 
ing the entropy constant fey, which is usually identified 
simply with Boltzmann's kg in the literature: the only 
certified fact one can be sure of is "fc^ — > kg for q — > 1" 
||, which entails that there is room to choose k? — k(q) 
in a suitable way. It is seen |l| that if one chooses (with 



Z q the generalized partition function) 



k T = k B z«-\ 



(2) 



k T q-l 



(1) 



in conjunction with the OLM formalism [0, the classical 
harmonic oscillator determines a specific heat C q = ks, 
so that the ordinary Boltzmann-Gibbs result arises with- 
out invoking the limit q — > 1. In addition, the OLM 
treatment is able to reproduce the zero-th law of ther- 
modynamics, a goal that had previously eluded Tsallis- 
thermostatistics practitioners p9| . One is then tempted 
to conjecture that many other ordinary results can be 
reproduced by Tsallis thermostatistics without invoking 
the limit q — > 1. In this paper, we revisit the ideal gas 
problem in Tsallis thermostatistics with such a goal in 
mind. More precisely, we perform an OLM treatment 
of the problem following the canonical ensemble struc- 
ture. We will show that the ordinary expression for the 
mean value of the energy is reproduced by Tsallis ther- 
mostatistics without the need of going to the limit q — > 1 , 
in contrast to the previous treatment of this problem. 

It is worth noting here that, as shown in Ref. p0[ , 
the usual thermodynamic limit N — > oo implies the limit 
q — > 1. There, the Tsallis probability distribution is de- 
duced in a way that mimics Gibbs' celebrated derivation 
of the canonical distribution for a system in contact with 
a heat bath [G0|. It was shown that Tsallis' distribu- 



1 



tion naturally arises for finite heat baths, the nonexten- 
sivity index q being related to the particle number TV 
that characterizes the bath. Gibbs' canonical distribu- 
tion, instead, results for infinite heat baths PQ]. This 
leads one to conclude that q goes over to unity in the 
limit iV — > oo. This result was found employing the un- 
normalized expectation values of Curado and Tsallis (see 
Ref. H ) . Analogously, if one uses the normalized expec- 
tation values of Ref. H , a similar argument leads to 



DN 



(I = 



DN -2 



(3) 



for the ideal gas in D dimensions. Since the limit N — > oo 
corresponds to q — ► 1, the same ordinary result is ob- 
tained from the nonextensive formalisms with various 
definitions of generalized expectation value (both nor- 
malized and unnormalized) in such a limit. 



II. BRIEF REVIEW OLM FORMALISM 

A general classical treatment requires consideration of 
the probability density /(x) that maximizes Tsallis' en- 
tropy, subject to the foreknowledge of the generalized 
expectation values of certain physical quantities. 

Tsallis' probability distribution [Q is obtained by fol- 
lowing the well known MaxEnt principle (^l). The 
Tsallis-Mendes-Plastino variational treatment Q in- 
volves a set of Lagrange multipliers Xj. The OLM 
technique developed in Ref. jl| pursues an alternative 
path that involves a different set of Lagrange multipli- 
ers, say, Aj : one maximizes Tsallis entropy S q in Eq. (|l|) 
^3| , ^Jl8| , subject to the modified constraints ("centered" 
generalized expectation values) |lp^|: 



I dx /(x) -1 = 0, 
/(x)« (o,(x) - ((0,}) 9 ) =0, 



(4) 

dx flxi'' (0;(xi ■ ■ "(),,,) = n. (5) 



where Oj(x) (j — 1,2, ...,M) denote the M relevant 
dynamical quantity (the observation level p2[). In the 
above, ((Oj)) are defined by [Q 



((Oj)) Q = 



jd* FtoQjW 

fdx/*(x) ■ 



(6) 



which are assumed to be a priori known. (The proce- 
dure given in Ref. (J] employs non-centered expectation 
values.) The resulting probability distribution reads |l| 



/(x) = Z- 1 



M 

l-ll-flE^O^-^)), 



(7) 



dx 



M 



l- { l-q)J2x' j (O j ( X )-((O j )) q 



(8) 



Although the procedure originally devised in Ref. j|] 
overcomes some problems posed by the old, unnormalized 
way of evaluating Tsallis' generalized expectation values 
^J,^3|, it yields probability distributions that are self- 
referential. The resulting distribution includes the inte- 
gral of the qth power of the distribution itself. This fact 
entails difficulties in numerical model calculations, for ex- 
ample. The complementary OLM treatment of Ref. fjj 
surmounts such difficulties. The above-mentioned self- 
reference problem does not arise in Eq. ((?]). 

It is shown in Ref. ]jj that the Lagrange multipliers 
Xj of the Tsallis-Mendes-Plastino procedure B and the 
corresponding A'- in OLM are connected to each other as 
follows: 



A' 



A; 



/dx /«(x)' 



(9) 



However, the genuine Lagrange multipliers are X'j in 
OLM. 

The probability density appearing in Eq. (0) is the one 
that maximizes the entropy S q which can be expressed 
in the alternative form fllfl 



zr 1 - 1 



q-l 

Also, the identical relation 

J dx /«(x) = Z\-« 
holds, from which it follows Q] that 

Sq = kx UXqZq, 

where ln x = (1 — x 1 ~ q )/(q — 1). 
Eq. (O) allows us to rewrite Eq. 



dx /»(x). 



as 



A' 



A, 



yl-q ' 



(10) 



(11) 



(12) 



(13) 



However, we again stress that the basic variables in OLM 
are A^-. 

Now, following Ref. Q, we define 

hXqZ' q =hXqZq-Y,^ ((0 3 )) q - 



(14) 



Introducing 



k' = k T z 1 -" 



(15) 



where Z q stands for the generalized partition function as in Ref. [Q, we straightforwardly obtain 



2 



dS q 
9{(O j )) q 



fc'A 



— (\n q Z>) =-<(0,-» ? 



(16) 
(17) 



Eqs. ([l6]) and (0) constitute the basic information- 
theoretic relations, on which statistical mechanics can be 
built a la Jaynes Q|. Here, one should remind the well 
known fact that in reconstructing statistical mechanics 
based on information theory, Jaynes could remove the 
concept of ensemble 21 1. In this sense, it is appropriate 
to emphasize that both the OLM jl| and Tsallis-Mendes- 
Plastino [Q formalisms employ identical a priori informa- 
tion, so that they are physically equivalent. 

Notice that k', as defined by Eq. (|l5|), obeys the con- 
dition k' — y &b as q — ► 1. It is a condition that this 
constant must necessarily fulfill. (See Ref. Jll|). Notice 
that if one makes use of the possib ility of choosing fcy as 
in one obtains, on account of Jig) , 



k' = k B , 



(18) 



i.e., the classical results are obtained without going to 
the limit q — > 1, as stated in the introduction. 

Moreover, another interesting result obtained from 
OLM is § 



fc'A^ = k T \ 



N 



2m' 



(23) 



where m is the particle mass, N the particle number and 
Pi the momentum of the ith particle. We are writing the 
TV-particle momenta collectively as P = (pi, P2, ■ ■ ■ , Pn)- 
One extremizes the entropy in Eq. (|l]), subject to the 
constraints [Q 



j dfi/(P) = 1, 

- c j dn/(p)«zr(p) = u< 



(24) 
(25) 



where dfl = (l/(N\h DN )) ]jf =1 dq t dp u with h the linear 
dimension (i.e. the size) of the elementary cell in phase 
space, and J Yii=i d<k = V N with the spatial volume V. 
We have also introduced the quantity c 



dfi/(P)«. (26) 



The equilibrium probability density is found to be 
f(P) = i- [1 - (1 _ q)((3/c)(H(P) - U q )}^ , (27) 

Z n 



J' 



(19) where the generalized partition function is 



which entails that the intensive variables are the same in 
both the OLM and Tsallis-Mendes-Plastino formalisms. 

As a special instance of Eqs. ([l6|), ( |l7| ) and (|l^), let us 
discuss the canonical ensemble. In this case, only a single 
constraint regarding the system Hamiltonian is consid- 
ered. Writing the associated Lagrange multiplier as (3' 
and the generalized internal energy as U q , 



where 



-^(ln q Z' g ) =-U q , 



In, Z' q = In, Z q - p'U t 



(20) 
(21) 

(22) 



The temperature T in Eq. (|20| ) is the same as that in 
the Tsallis-Mendes-Plastino formalism. 



III. IDEAL GAS IN 
TSALLIS-MENDES-PLASTINO FORMALISM 

The classical ideal gas in D-dimensional space in the 
Tsallis-Mendes-Plastino formalism has been considered 
in Refs. (^||. For comparison, we recapitulate here some 
of its main results in the case < q < 1. The Hamilto- 
nian reads 



Z q = I <ffi [1 - (1 - q)(J3/c)(H(P) - u q )\ — 



and we have 



c = [Z q ] 



(28) 



(29) 



In the above, (3 is the Lagrange multiplier associ- 
ated with the constraint in Eq. (25). Notice that Eq. 
( p7| ) is indeed a self-referential expression through c. In 
the present circumstances, however, this problem can be 
overcome, since a straightforward mathematical manip- 
ulation yields 



pu q 



DN 



(30) 



with c given by 



1-9 



V 



^1 + dn_} Nlh DN 



l + (l-g) 



DN 



2irm 



2(l-q) 

q(l — q)+DN/2] l-(l-«)b» 



(31) 



To discuss the statistical properties of the particle en- 
ergies, the ith and jth single-particle Hamiltonians may 
be considered, namely, Hi = pi /2m and Hj = p|/2m, 



3 



respectively. Their generalized variance, covariance and 
correlation coefficient are defined by || 

(A g ^) 2 = ((^ 2 )} g -((^}}2, 
C^HuHi) = ((HiHj)), - ((H^dHj)),, (32) 

C q (Hi,Hj) 



p(H tl H 3 



V(A 9 ff 4 ) 2 (A^) 2 



respectively. A straightforward calculation shows that 
these quantities are given by 



(AgHi) = 



c 2 2D+(l-q)D 2 {N - 1) 



2(3 2 4 - 2q - (1 - q)DN ' 



C q {Hi, Hj) 
p{H l ,H ) = - 



(l-q)D 2 



2(3 2 4-2g+ (1 - g)D7V' 
(l-q)D 



(33) 



2 + (1 - q)D(N - 1) 
In the limit q — * 1, one gets 

C q {H h Hj)^Q, 



(34) 



IV. IDEAL GAS IN OLM FORMALISM 

We revisit here the classical ideal gas problem consid- 
ered in the previous section within the OLM framework. 
We extremize Eq. (|l]), subject now to the constraints 



J dflf(P) -1 = 
J dQ/(P) 9 (ff(P) - U v ) = 0, 



(35) 
(36) 



with H given in Eq. (|23|). 

The probability distribution /(P) reads 



/(P) = ^[1 - (1 - g)/3'(ff(P) - (7 g )] — , 



where 



= / dfi[l - (1 - q)p'(H(P) - U q ) 



(37) 



(38) 



In these equations, (3 1 is the La gran ge multiplier asso- 
ciated with the constraint in Eq. (|36|). 

We now follow the steps indicated in Ref. ||] and sum- 
marized in Section [ll| Our interest lies in the interval 
< q < 1 again. First, define 



Ri = 



p| 2-9 

' 1-9 



/ 2-g , DN_ 
\l-q 2 



(39) 



and 

R-2 = 



27rmc 


2 


[(1 -?)/?'. 





l + (l-q) 



(3'U q 



(40) 



The associated generalized partition function is given 



by 

Z q {p') = Rx 
while, introducing 

Gi 



V 



N 



Ri 



DN 



20>Z q (p>y 



the generalized internal energy turns out to be 

V N 



Uq — G\ R\ 



R 2 . 



Nlh DN 

After replacing Z q ((3') by Eq. (gl|), one finds 

_DN 
Uq ~ 2/3' ■ 



(41) 



(42) 



(43) 



(44) 



Now, taking advantage of the fact that (Cf. Eq. (|2' 

1 



k R T' 



one is immediately led to 



Un = ^Nk B T, 



(45) 



(46) 



i.e., to the ordinary result. Note that this result is ob- 
tained here without going to the limit q — > 1. By following 
the methodology advanced in Ref. we see that, with 
the identification of Eq. (|^), Tsallis thermostatistics is 
able to reproduce ordinary results independently of the 
specific choice of the value of q. 

Now, let us discuss statistical properties of the single- 
particle energies. Following again Section ||, we focus 
our attention on the single particle Hamiltonians per- 
taining to the ith and the jth particles, Hi = p 2 /2m and 
Hj = p 2 /2m, respectively. Their generalized variance, 
covariance and correlation coefficient are defined collec- 
tively in Eq. (yields 



(A^) 2 
C q (Hi, Hj) 



1 2D + (1- q)D 2 {N -1) 
2(^') 2 A-2q-{l-q)DN ' 

1 (l~g)£ 2 

2{f3') 2 A-2q- {\- q)DN' 
(l-q)D 



(47) 



p{H u Hj)- 2 + { i_ q)D ^ N _ x y 
respectively. In the limit q — > 1, one has 



4 



C q (Hi, Hj)—* 0, 



(48) 



We see that the limit g — ► 1, C g and /? behave as in 
the Tsallis-Mendes-Plastino formalism in Ref. ||. How- 
ever, for (A q Hi) 2 , we find something new: the result be- 
comes independent of the density. It depends only on the 
temperature, as it does in the ordinary Boltzmann-Gibbs 
treatment. 



V. CONCLUSIONS 

In this work, we have addressed the ideal gas problem 
within the framework of Tsallis thermostatistics reformu- 
lated by the optimal Lagrange multipliers (OLM) method 
advanced in Ref. Q . Three interesting observations have 
been presented: 



• The internal energy was found to be given by 
U q = ^Nk B T, 



which is the same as the ordinary result obtained 
by Boltzmann-Gibbs theory. We derived this result 
without going to the limit q — > 1. That is, it is valid 
for all q. 

• The correlation coefficient for the particle energies 
become density-independent in OLM, in contrast 
to that in the Tsallis-Mendes-Plastino formalism 
in Ref. J§. 

• The constant fcr in the definition of the Tsallis en- 
tropy may depend on q and, therefore, on the num- 
ber of particles. In the Boltzmann-Gibbs canoni- 
cal ensemble theory, the limit particle number N 
going to infinity is always in mind, both explic- 
itly and implicitly. This indicates that kx — > ks 
and q — ► 1 may correspond to N — ► oo w 
Avogadro' s number. 
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